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Spatial 't Hooft loops of strength k measure the qualitative change in the behaviour 
of electric colour flux in confined and deconfined phase of SU(N) gauge theory. 
They show an area law in the deconfined phase, known analytically to two loop 
order with a "k-scaling" law k(N — k). In this paper we compute the 0(g 3 ) 
correction to the tension. It is due to neutral gluon fields that get their mass 
through interaction with the wall. The simple k-scaling is lost in cubic order. The 
generic problem of non-convexity shows up in this order. The result for large N is 
explicitely given. 



1. Introduction 

The spatial 't Hooft loop 1 monitors the behaviour of chromoelectric flux. 
In the confined phase of QCD the flux behaves quite different from that 
in the deconfined phase. Hence it is an order parameter for gluodynamics, 
and to some extent also for full QCD 3 . 

It is a perturbatively calculable quantity . More precisely it can be 
computed as a tunneling effect through a potential mountain that is per- 
turbatively calculable. The potential has a periodic structure due to Z(N) 
symmetry 1 . This Z(N) symmetry shows up in every gauge system with a 
compact periodic dimension and leads to an area law for the loop at high 
temperature. Hence the calculation described below is equally valid in for 
domainwalls one will encounter in a four dimensional world with a compact 
fifth dimension . 

In this note we compute the 0(g 3 ) correction to the tension in the large 
TV limit. 
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2. What is the 't Hooft loop? 

The 't Hooft loop Vk(L) is usually defined 1 as a closed magnetic flux loop 
L of strength fc^, with N the number of colours. In operator language 
it is defined as a gauge transformation with a discontinuity expik^- when 
crossing the minimal surface. As we will justify below, one can write this 
gauge transform for L in the x-y plane as a dipole sheet 3 in the following 
way: 

47T f 

V k {L) = expi— / dxdyTrE z Y k . (1) 

iV JS(L) 

The NxN diagonal traceless matrix Y& is defined as 

y fe = diag(k, k...k, k-N,k- N , k — N) (2) 

with N-k entries k and k entries k-N, to have a traceless matrix. The charges 
Yfe are generalizations of the familiar hypercharge, with k=l. The charge 
Yfe of a gluon is or ±N. The multiplicity of the value N is k(N — k). The 
same is true for the value — N. So, e.g. for N = 3 and k = 1 one finds the 
four kaons with hypercharge ±3. 

Exponentiation of Yfe gives expi^Y^ = cxpikj^- = Zk, the centergroup 
element. 

E z is the z component of the canonical electric field strength operator 
E = \ a E a . 

3. How to compute its thermal average. 

Apart from its obvious connection with colour electric flux the 't Hooft 
loop is intimately related to the Z(N) symmetry. Take its order parameter 
P(Aq), and move it through the dipole layer. Then it will get multiplied 
by the discontinuity expik^-, being a fundamental test charge. Immersing 
the loop in the plasma induces a disturbance. The disturbance is described 
by a profile C. Since the loop is gauge invariant the response of the plasma 
is too. This profile is the phase of the Polyakov loop as a function of its 
distance to the minimal surface. 

So the approach will be to compute the free energy excess AF due to 
the presence of the Polyakov loop profile. Let the box be of size L\ r * L z , 
with L z >> L tr and both macroscopic. Extend the loop to the full x-y 
cross-section of size L\ r and located at say z = 0. 



a The A matrices being normalized to Tr \ a \b = 2^ a . 6 ' 
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Then we have exp-AP(C) = / DCcxp -^fU(C). U{C) is the con- 
strained effective potentiai: 

exp-^£/(P) = J DA DAS(P - P(A )) exp ^S(A,s k ). (3) 

and we have used an abbreviated notation for the constraint: 

5(P - P(Ao)) = U Itl 5(pW(z) - P l (A a (z))) (4) 

where I runs from 1 to N — 1. In the following we parametrize the fixed 
loop by a diagonal traceless NxN matrix C(z) = diag(Ci,C2, Cjv) : 

p(i) = —TrexpilC(z). (5) 

So the effective potential is defined on the Cartan space in which the ma- 
trices C live. The matrix describes the profile of the loop. 

It also fixes the path in Cartan space (the space of diagonal C's) along 
which the minimal profile is realized. The path turns out to be the sim- 
plest possible 5 : if parametrized by q, < q < 1, it is given by the one 
dimensional set of Cartan matrices 

Y k (q) = qY k (6) 

with Yfc the charge characterizing the strength of the dipole layer, eq.(l). 
In exponentiated form it goes from 1 to expifc^, as q goes from to 1. 
Z(N) invariance of the profile functional U(C) garantuees we we can take 
a smooth path from 1 to expi/c^, instead of a path that makes a jump 
expiA;^ at the loop and returns to 1. A proof of the rectilinear path being 
the minimal one is still lacking. Only for SU(3) and SU(4) it is known to 
be the case by inspection 4 5 and a proof at large N is given in ref. 4 . 

4. Perturbative results 

A loop expansion of the constrained potential eq. (3) is straightforward. One 
introduces the colour diagonal background field B through A^ = P<5 M ,o + 

The saddle point expension around the Polyakov loop P(C) fixes the 
background B in terms of C . 

The results from one and two loops are known 4 5 and are given by 

4-7T 2 

Pk (L) = k(N k) 2M ^ T\1 1.1112.. (a s N) + 0(g 3 )) (7) 
The gauge coupling is running according to the MS scheme 8 . 
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5. Three loop potential and self energy matrix 




Figure 1. The only three loop diagram of free energy topology with an infrared diver- 
gence. The shaded blob is the one loop selfenergy of. The colour index c, d of the two 
propagators need not be the same due to background dependence inside the blobs. 

This diagram is in the case of vanishing background linearly divergent 
for the Coulombic components of the propagator. For notational conve- 
nience we write II^, (q) for the Coulombic selfenergy matrix at zero mo- 
mentum and frequency, dropping the Lorentz indices. The index k means 
we are evaluating the selfenergy along the tunneling path Y}-{q). It is not 
diagonal except for k = 1. The eigenvalues Il{ k \q) 2 turn out to be of the 
form: 

l-6^(q-q 2 ) = (f(r,q) + l) i 
So the inverse Coulomb propagator in the background q becomes: 

«Xc + ng(?) (8) 

and gives an order g 3 contribution to the effective potential. In particular 
for q = 0, this contribution is proportional to m 3 D . 

For the infrared finite result to be Z(N) invariant we should keep con- 
tributions from all Matsubara frequencies in the selfenergy. 

Subtract from from all sclfcncrgics in 3 or more loops the zero frequency- 
momentum piece, to get convergent integarls. Divergencies due to the mag- 
netic sector remain. 

6. The 0(g 3 ) contribution to the tension. 

We have now eliminated the infrared divergencies from the potential up and 
including three loop order. By rcsumming the Z(N) invariant eigenvalues 
this procedure respects the Z(N) invariance of the potential. 
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The price to pay is a negative selfenergy in a certain window of q values 
of order 1. The problem has not to do with the infrared scales. It is a 
generic problem occurring in quantum corrections to tunneling through a 
barrier. However, for a window of loop strengths i < r = -j| < | the mass 
stays positive over the whole range of q values. In particular in the large N 
limit we will be able to extract the potential without this problem. 

The cubic term in the tension is obtained by minimizing the effective 
potential cq.(3). 

The tension becomes then to cubic order: 

p r (T) = p« (1 - 1.1112..(a s AT) + ( A)l/ (r)(asA r)f + (a 2 s )) (9) 

and 

w J 6q(q - l)r(l - r) 1 ; 

As in the pressure 8 it contributes with a sign opposite to the 0(g 2 ) term. 

The attraction between loops becomes stronger due to the convexity of 
I(r). To see this, take the ratio of the tension of the k-loop and compare 
it to the k times the tension of the elementary loop with k=l: 

PJl(T) 3 , i 

= (1 - r)(l + (^)Hl(r) J(0))(a 8 JV)* + 0(a 2 ) + O(^) (11) 

The cubic correction has now a negative coefficient, so that the ratio is 
smaller due to the presence of this correction. 
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